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In this talk we review the status concerning vacuum integrals needed in perturbative expansions of QCD at 
non-zero temperature. We will focus on the differences as compared to familiar zero-temperature techniques, and 
provide a list of known basic master integrals. 
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1. Introduction 

Weak-coupling expansions in finite tempera- 
ture field theory can be organized in a way 
that is very similar to the better-known zero- 
temperature case. In particular, for observables 
describing the physics of thermally equilibrated 
systems, the so-called imaginary time formalism 
provides a very close analogy. 

For a given system that is coupled to an exter- 
nal heat-bath with temperature T, the relevant 
observables have to be averaged with the density 
matrix, (O) = Tt(Oc- h / t ) /Tr(e~ H ^ T ), where 
H denotes the Hamiltonian of the system. In the 
case of interacting quantum fields <J) in thermal 
equilibrium, one can lift this to a functional inte- 
gral representation of the averages, 
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with Euclidean action and where the trace is en- 
forced by constraining the path integral to fields 
that are periodic in imaginary time r: 



+ 1/T,x) =±<p{r,x) 



(2) 



where the upper/lower sign refers to 
bosonic/fcrmionic fields. 

Finite-temperature Feynman rules hence dif- 
fer from zero-temperature ones by the fact that, 
due to the compact support in imaginary time r, 
the integration measure includes a discrete sum 
J -&- — > T y\ over the Matsubara frequencies 
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Po = litTn for bosons and po = 2irT{n + i) 
fcrmions, respectively, where bEZ. As the pres- 
ence of a preferred direction - the rest frame of 
the heat-bath - has broken the 4d Lorentz sym- 
metry, one works with four-momenta P = (po , p) 
and uses Euclidean metric g^ v = 5^, such that 

P 2 =pI+p 2 - 

In this short note, we will assume the reader 
be familiar with standard methods of zero- 
temperature perturbation theory, such as inte- 
gration by parts (IBP) p] and its systematic al- 
gorithmic use [2], and then point out the main 
differences that occur in a thermal system. For 
illustration purposes, we will concentrate on the 
most basic class of sum-integrals, vacuum (or 
bubble) integrals, which are needed to compute 
e.g. thermodynamic observables that directly re- 
late to vacuum diagrams, like the pressure (see 
3 4 and references therein). These vacuum in- 
tegrals are at the same time needed as building 
blocks for higher-point Greens functions, after ex- 
panding in external momenta, see e.g. [5]. A col- 
lection of analytically known integrals that are 
relevant in these contexts will be presented. 



2. Bosonic master integrals 

We work in dimensional regularization, and for 
convenience take the integral measure to be 
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2.1. One-loop bosonic integrals 

The basic bosonic 1-loop tadpole integral is 
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(P 2 )™ 
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(27rT) 2 «- m r(i)r(n) 
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It obeys the recursion relation Z^!S = 
2 "~2n +2e -^™' as can either be derived from in- 
tegration by parts (IBP) in the integral over 3- 
momentum, or confirmed directly from the above 
explicit solution. Note that, in contrast to the 
zero-temperature case, already at 1-loop level one 
finds in principle infinitely many master integrals. 
In practice, however, only a small number of them 
appear in any practical computation. 

1-loop master integrals that are needed for the 
pressure are 

T 2 4 e r(l + £ ) C(~l + 2e) 1 
12 r(i) C(-l) l-2e 

« ^[l + e(2- 7E + 2Z 1 ) + ...] , 
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where Zj = ^pjp 
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2.2. Two-loop bosonic integrals 

At 2-loop order, there is a major simplification: 
the basic sunset-type master vacuum integral (cf. 
Fig. [T]) vanishes identically, 
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as can be proven by IBP. In fact, using that the 
integral of a total derivative vanishes in dimen- 
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Figure 1. 2-, 3- and J^-loop bosonic master sum- 
integrals. Each solid line corresponds to a bosonic 
propagator l/P 2 , each dashed line to an inverse 
propagator P 2 . A dot on a line denotes the same 
object squared. 



sional regularization 
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and choosing (here, d = 4 — 2e denotes the space- 
time dimension) 



+ 752 (P° + 1o)(poqt - qoPi) , 
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fi = (d-3)(pj 

2 

Q 2 " 

after working out the derivatives, using the shift 
Q —> P — Q, and exploiting symmetry of the sum- 
integral under P <-^ Q, one gets 

= (d-3)(d-4)5, 

which completes the proof of Eq. ([5]). 

2.3. Three-loop bosonic integrals 

At the 3-loop level, a number of master inte- 
grals have been encountered in the computation 
of the pressure. For a classification, let us adopt 
a naming scheme advocated in [3J, 
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Note that M. is symmetric in its two indices. The 
known cases (cf. Fig. [I]) include a basketball-type 
and a spectacles-type integral 

Mo ° " (4^F 24^[ 1 + e6ll + £26l2 + -] ' (9) 
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as well as the modified basketball-type integral 
that is needed at higher orders [I], 

M' = U - 
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The expansion coefficients read |6I4) 

bu = f§ - 3 7e + 8Zi - 2Z 3 , (12) 

621 - f+ 7E + 2Z!, (13) 

622 = 48.7976.. , (14) 

"Hi = I - TT7b + ff ^1 - 7^3 , (15) 

while 612 has not been computed yet. 

Other integrals of the class Eq. ([8]) can be re- 
duced by standard methods, like shifts of inte- 
gration momenta exploiting the symmetry of the 
integrand, such as 

M-u = -^Moq + 22? 5, (16) 

or, in a graphical representation 

where the notation is like in Fig. [TJ 

2.4. Four-loop bosonic integrals 

At the 4-loop level, the classification and eval- 
uation of sum-integrals has not been tackled in a 
systematic way yet. There is however some pio- 
neering work exploring scalar theory at this or- 
der [4]. The major simplification (as compared 
to QCD) is that the reduction step is trivial, and 
there is only one genuine 4-loop sum-integral to 
compute. Hence, the only 4-loop master integral 
that is known presently is (cf. Fig. [I} 
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with coefficients 
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tn = f - 4 7e + 12Zi - 4Z 3 - 3C(2) , (18) 
t 12 = 2612 + 25.7055.. - 3C(2)(28.9250..) . (19) 



These coefficients have been calculated in [3], 
profiting from the particular structure of the in- 
tegral Eq. (|17|) . which contains three insertions 
of 1-loop bubble-type integrals, and can hence be 
tackled in complete analogy to the one in Eq. ([9]) , 
as pioneered in [6]. The numerical values given 
above result from simple numerical integrations 
and are known to more than 10 digit accuracy. 

Note that in the physics calculation where this 
integral was needed, it occurred in combination 
with the one in Eq. ([9]) such that the number 612 
cancels in the final result. 

3. Fermionic master integrals 

Let us denote fermionic four-momenta by 
braces, {P}, to indicate that their po = 2-kT(ji + 
i) with n G Z. 

There is an important class of relations between 
fermionic and bosonic sum-integrals, which can 
be derived by partitioning the Matsubara sums 

E=E+E ( 2 o) 

ii£Z n even n odd 

and then rescaling spatial integration momenta 
on the left-hand-side as pi — > \pi (see, e.g., [7]). 

3.1. One- and two-loop fermionic integrals 

The relations obtained as just described reduce 
all fermionic 1- and 2-loop tadpoles to the respec- 
tive bosonic case, 
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where in the last step Eq. ([5]) was used. 



3.2. Three-loop fermionic integrals 

To classify the master integrals at the 3-loop 
level, let us again adopt the notation of [3], 
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Figure 2. 2- and 3-loop fermionic master sum- 
integrals. Each arrow-line corresponds to a 
fermionic propagator. 



Figure 3. Graphical representation of a linear re- 
lation between various J^-loop sum-integrals. No- 
tation is as above. 
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Note that, while M is still symmetric in its two 
indices, M is not. 

It turns out that, by the same partitioning 
trick as above, the two possible fermionic 3- 
loop basketball-type vacuum integrals are related, 
such that it suffices to include one of them into 
the basis of master integrals, 
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where [5] 

&31 = W- 3 7b-# ln2 + 8Zi-2Z 3 . (27) 
Further 3-loop masters (see also Fig. [2]) arc 
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In close analogy to the bosonic case Eq. fT6J), 
others get reduced by applying a shift of integra- 
tion variables and exploiting the integral's sym- 
metry, 



A/Lu 
M-n 



-\n 00 + 21? 5, 
-iMoo+Z?S+??S 



(34) 
(35) 



3.3. Four-loop fermionic integrals 

Not a single genuine fermionic 4-loop sum- 
integral has been computed to date. However, 
at this level, a theoretically important contribu- 
tion is expected to be made for the QCD pressure 
[§], adding the last missing perturbative piece of 
a physically complete leading-order description, 
which includes (known) non-perturbative coeffi- 
cients. Since, naturally, an analytic treatment of 
diagrams containing quark loops seems to be a 
good starting point (due to the simpler structure 
of fermionic vertices and propagators), progress 
with this type of sum-integrals can be expected 
in the near future, before the gluonic sector (cf. 
Sec. 12. 4p will be classified. 

To obtain a first useful relation, one can apply 
the strategy of Eq. (j2"0|) to the bosonic integral of 
Eq. (flT)) . to get a linear relation between the four 
distinct fermionic integrals of the same topology, 
as depicted in Fig. [3J 

4. Chemical potential 

In systems which allow for an excess of, say, 
particles over anti-particles, one introduces fur- 
ther parameters - the chemical potentials fi - as 
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Lagrange multipliers of the particle number oper- 
ators into the Lagrangian. The effect is to change 
four-momentum to be used in the free propaga- 
tors of the respective particle species as 



Po 



Po -l£t 



(36) 



such that all sum-integrals pick up a dependence 
on the parameter(s) /i. In the case of QCD, to 
accommodate for a non-zero baryon number den- 
sity, chemical potentials have to be introduced for 
all flavors of quarks. 

In fact, contributions to the QCD pressure have 
been calculated including the effects of a quark 
chemical potential, requiring a generalization of 
the sum-integrals discussed in Sec. [3] to include 
the /i-dependence. This was done through 3- 
loop order in [10j . where all required expansions 
are shown. While for the 1-loop fermionic tad- 
pole an analytic result can simply be given in 
terms of the generalized Zeta function £(x, z) (cf. 
Eq. (|40|) below), the 3-loop generalizations A/y (/x) 
and M-ij(n) are somewhat more complicated to 
evaluate, and their expansions turn out to contain 
derivatives d x C,{x, z) [TU] , 

One notable structural change as compared to 
the case /j, = is that certain relations that relied 
on symmetries of the integrand no longer hold 
true. In particular, there are now non- vanishing 
/i-dependent 2-loop sum-integrals of sunset-type, 
such as (cf. Eq. @) 
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where z = | — i^pf an< ^ C( x i z ) — 9 x ((x,z). 
Note that the result Eq. ([37]) is compatible with 
Eq. (|22|) in the limit /i — ► 0, as it has to be. 

We will not list the further known integrals 
here, but finally briefly discuss a different setting, 
in which related sum-integrals appear. 



5. g-integrals 

In some cases, one might be interested in in- 
tegrals which (for fj, = 0) interpolate smoothly 
between the bosonic and fermionic cases, using 

Po = 2TrT(n + q) , n e Z and q e R. . (39) 

Below, four-momenta with such q-dependent 
components will be denoted by brackets, [P]. 
Hence, all sum-integrals pick up a dependence 
on the parameter q, and, since the integers n 
are summed, are periodic with respect to a unit 
change in q, such that it is sufficient to determine 
them on the interval q G [0, 1). Note that q = 
corresponds to the bosonic integrals discussed in 
Sec. [21 while q — i corresponds to the fermionic 
ones discussed in Sec. [3] 

Notably, the need for such g-dependent sum- 
integrals arises in attempts to analyze the be- 
havior of the so-called spatial 't Hooft loop [12] , 
which can be taken as an order parameter for 
the deconfincment phase transition in pure Yang- 
Mills gauge theory, and requires the minimization 
of a certain effective action with respect to the 
parameter q |13j . 

To show a concrete example, the basic bosonic 
1-loop tadpole can now be represented as 
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x [(-l) m C(2n-m-3 + 2e,l-g) + 
+ C(2n - m - 3 + 2e, q)} , (40) 

where ({x,q) = E T T=o («+V is thc g eneral " 
ized Zeta function, with £(ir, 1) = £(x) and 
C(x, i) = (2 X - l)C(x) relating J™(g) to Eq. Q 
and Eq. pT|) . 

We will not discuss higher loop orders here, al- 
though some results are also available in the liter- 
ature [13] , while others (notably up to the 3-loop 
level) can in principle be deduced from [10111] by 
analytically continuing the g-dependence to com- 
plex values q — > z = \ — in/(2nT), such as the 
g-dependent 2-loop sunset-type sum-integral 

T 2 (2q-l) 2 
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as obtained by explicit calculation, or by analytic 
continuation from Eq. ((37|) . It vanishes as q — * | , 
as it should (cf. Eq. (|2"2"]1). 
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